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Abstract: 

Prion diseases, traditionally referred to as transmissible spongiform encephalopathies 
(TSEs), are invariably fatal and highly infectious neurodegenerative diseases that affect 
a wide variety of mammalian species, manifesting as scrapie in sheep and goats, bovine 
spongiform encephalopathy (or mad-cow disease) in cattle, chronic wasting disease in 
deer and elk, and Creutzfeldt-Jokob diseases, Gerstmann-Strussler-Scheinker syndrome, 
W I fatal familial insomnia, and kulu in humans, etc. These neurodegenerative diseases are 

O ' caused by the conversion from a soluble normal cellular prion protein (PrP*^) into insol- 

r^ ■ uble abnormally folded infectious prions (PrP^'^) and the conversion of PrP*-^ to PrP^'^ is 

believed to involve conformational change from a predominantly a-helical protein to one 
rich in /3-sheet structure. Such conformational changes may be amenable to study by 
molecular dynamics (MD) techniques. For rabbits, classical studies show they have a low 
^ . susceptibility to be infected by PrP^'^, but recently it was reported that rabbit prion can 

^^ I be generated through PMCA in vitro and the rabbit prion is infectious and transmissi- 

2 ■ ble (Proc Natl Acad Sci USA 109(13) 5080-5085 (2012)). This paper studies the NMR 

f"^ . structure of rabbit prion protein wild-type and mutants by MD techniques, in order to un- 

T^J- ' derstand the specific mechanism of rabbit prion protein and rabbit prions. 

o 

Keywords: prion diseases; rabbit prion proteins; rabbit prions; NMR structures; molecu- 
lar dynamics study 



1 Introduction 

"Mad cow disease" (Bovine Spongiform Encephalopathy (BSE)) belongs to a contagious 
type of Transmissible Spongiform Encephalopathy (TSE). Scientists believe it is caused 
by Prions (the misfolding prion proteins) but they may have not yet solved the riddle 
of "mad cow disease". This is due to a Prion is neither a virus, a bacteria nor any mi- 
croorganism so the disease cannot be caused by the vigilance of the organism immune 
system and it can freely spread from one species to another species. The humans exists 
the susceptibility of TSEs. For example, the human version of mad cow disease named 



Creutzfeldt-Jakob Disease (CJD) and variant CJD (vCJD) just happen randomly through 
infections of transplanted tissue or blood transfusions or consumption of infected beef 
products. Cats, mink, deer, elk, sheep and many other animals are also susceptible to 
TSEs. However, rabbits, horses and dogs seem to be unaffected by Prions lfT6l 171 [T3l. 
Scientists do not know the reason. 

The prion protein is a naturally occurring protein in vivo. Its lesions in brain are not 
caused by the vigilance of the immune system. Recent studies have found that the lesions 
led astray as long as by contact with other normal prion proteins. The cells are arranged in 
accordance with the instruction of the gene and formed into proteins with different shapes 
and functions. But, like cardboard boxes, proteins need to be properly "folded" in order 
to ensure their normal work. When proteins are folded into the wrong shape, they do 
not work. Under normal circumstances, the cells will supervise these misfolded proteins 
and automatically decompose them. However, the supervision mechanism is not with 
100% insurance. Scientists found that the rate of decomposition of Prions is not quick 
enough and these Prions accumulate and change the cellular metabolism and eventually 
kill the cell. This leads to the death of neurons in the brain. The dead neurons decompose 
and release more prion proteins into the biological mechanism to cause prion diseases. 
The infectious diseased prion is thought to be an abnormally folded isoform (PrP^'^) of a 
host protein known as the prion protein (PrP*-). The conversion of PrP*^ to PrP^'^ occurs 
post-translationally and involves conformational change from a predominantly a-helical 
protein to one rich in /3-sheet amyloid fibrils. Much remains to be understood about 
how the normal cellular isoform of the prion protein PrP*^ undergoes structural changes 
to become the disease associated amyloid fibril form PrP^'^. The "structural conforma- 
tional" changes of PrP^'^ from PrP*-^ are just very proper to be studied by MD techniques. 
Classical studies have showed that rabbits have low susceptibility to be infected by PrP^'^ 
IDIllIllSIBlIISllTliniinillllEEini- However, recently there is one exception saying 
that (i) rabbit prion can be produced through saPMCA (serial automated Protein Misfold- 
ing Cyclic Amplification) in vitro I2j^4| (though not by challenging rabbits directly in vivo 
with other known prion strains), and (ii) the rabbit prion generated is infectious and trans- 
missible |2|. Rabbit prion protein had already had NMR IfTOl and X-ray 171 structures 
in Protein Data Bank (www.rcsb.org) with PDB entries 2FJ3, 3079 respectively. This 
paper will study the NMR structure of rabbit prion protein (and its I214V and S 173N mu- 
tants: NMR structure PDB entries 2J0M, 2J0H) by MD techniques, in order to explain 
the specific mechanism about rabbit PrP*^ (RaPrP*-) and the conversion of PrP^ — )-PrP^'^ 
of rabbits. Surely, the studies should provide valuable knowledge about the rules gov- 
erning the PrP*-^ — >PrP^'^ conversion, which can provide some ideas for designing novel 
therapeutic approaches that block the conversion and disease propagation. 

The rest of this paper is organized as follows. In Section 2, we will briefly review 
previous research results on rabbit prion protein listed in H] |5] [8] [3] \T5l [161 El [El [111 
[T8l|2Tl[TTl|2l|4|. Section 3 will present the MD simulation materials and methods for 
NMR structures of RaPrP*^ wild-type and mutants. Section 4 will give MD simulation 
preliminary results and their preliminary discussions. At last, concluding remarks on 
RaPrP are summarized in Section 5. 



2 Brief Review on Rabbit PrP 

The symptoms of TSEs were first described for sheep in 1730 and called "scrapie" in 
England, "vertige" in France and "Traberkrankheit" in Germany (T4\ . Now we know 
that many species such as sheep, goats, mice, humans, chimpanzees, hamsters, cattle, 
elks, deers, minks, cats, chicken, pigs, turtles, etc are susceptible to TSEs. But many 
laboratory experiments show that rabbits, horses and dogs seem to be the resistant (or 
at least the low-susceptibility-rate) species to TSEs. However, recently Chianini et al. 
(2012) reported "rabbits are not resistant to prion infection" Q. Thus, at this moment it 
is very necessary for us to briefly review some laboratory works on RaPrP: 

• Barlow and Rennie (1976) made many attempts to infect rabbits with the ME7 
prion strain and other known prion strains but all failed at last HI. 

• Korth et al. (1997) found that RaPrP was not recognized by a conformational 
antibody specific for PrP^'^-like structures |j8]. 

• Vorberg et al. (2003) found multiple amino acid residues (such as GLY99, MET108, 
SER173, ILE214) within the rabbit prion protein inhibit formation of its abnormal 
isoform [i6|. The authors made some substitutions of mouse PrP amino acid se- 
quence by rabbit PrP amino acid sequence and found (i) at the N-terminal region 
(residues 1-111) the PrP^'^ formation is totally prevented, (ii) at the central region 
(residues 112-177), the constructed PrP failed to be converted to protease resis- 
tance, (iii) at the C-terminal region (residues 178-254) the formation of PrP^'^ is 
drastically decreased but is not abolished completely. Thus, rabbit cells are nega- 
tively affected by the formation of PrP^'^. 

• Femandez-Funez et al. (2009) found RaPrP does not induce neurodegeneration in 
the brains of transgenic flies ||5]. 

• Nisbet et al. (2010) created a mutant mouse PrP model containing rabbit PrP- 
specific amino acids at the GPI anchor site and found that the GPI anchor attach- 
ment site (w site) controls the ability of PrP*^ ^^PrP^'^ and the residues at uj and 
uj+l of PrP are important modulators of this pathogenic process llT2l . 

• Khan et al. (2010) found the propensity to form /3-state (the /3-sheet-rich structure) 
is greatest for hamster PrP, less for mouse PrP, but least for the PrP of rabbits, 
horses and dogs under different conditions and using two-wavelength CD (Circular 
Dichroism) method they also found a key hydrophobic staple-like helix-capping 
motif keeping the stability of RaPrP's X-ray crystallographic molecular structure 

m. 

• Wen et al. (2010) using multidimensional heteronuclear NMR techniques reported 
that the 12 14V and S173N substitutions result in distinct structural changes for 
RaPrP*-^ [17, 18|. The authors also reported the highly ordered /32-a2 loop may 
contribute to the local as well as global stability of the rabbit PrP protein lITSll . 



• Zhou et al. (2011) found that the crowded physiological agents FicoU 70 and 
dextran 70 have effects significantly inhibiting fibrillation of rabbit prion protein 

midii. 

• Chianini et al. (2012) generated rabbit PrP^'^ in vitro subjecting unseeded normal 
rabbit brain homogenate to serial automated protein misfolding cyclic amplification 
(saPMCA) and found the rabbit PrP^'^ generated in vitro is infectious and transmis- 
sible 12 1 . We noticed that the rabbit prion is just produced through PMC A in vitro 
not by challenging rabbits directly in vivo with other known prion strains. 

3 Materials and Methods 

This is the continuation and extension of the MD simulations of [19|. Zhang (2010) 
carried out (i) 15 ns of production phase of MD simulations and (ii) the heating phase 
of MD simulations is starting from 100 K with one set of initial velocity (denoted as 
seed2) II19L This paper continued to finish another 15 ns of production phase of MD 
simulations for seed2, and carried out other two sets of initial velocities for heating from 
100 K (denoted as seedl and seed3) of MD simulations with 30 ns of production phases. 

4 Preliminary Results and Preliminary Discussions 

For seedl'^seedS, during the whole 30 ns, radii of gyrations of the wild-type and the 
I214V and S173N mutants have been always level off around 15 angstroms whether in 
neutral pH environment or in low pH environment. 

In neutral pH environment, compared with the 12 14V and S173N mutants, the wild- 
type is always at the lowest level of RMSD (root mean square deviation), RMSF (root 
mean square fluctuation) and B-factor values. This show to us the wild-type RaPrP is 
very stable under neutral pH environment. However, under low pH environment, we 
cannot see some large differences among the wild-type, the 12 14V mutant and the S173N 
mutant. 

The variations of the molecular structures of the wild-type and its I214V and S173N 
mutants during the 30 ns for the three sets of MD simulations can be seen in Figures 
1~3, where the graphs were produced in the use of DSSP program |0, H is the a-helix, 
I is the TT-helix, G is the 3-helix or 3/10 helix, B is the residue in isolated /? -bridge, E 
is the extended strand (participates in /3-ladder), T is the hydrogen bonded turn, and S 
is the bend; HI, H2, and H3 respectively denote the a-helix 1, a-helix 2 and a-helix 3 
of a prion protein. We can see that for the wild-type the three a-helices have been still 
kept during the whole 30 ns of each set under the neutral pH environment, but under low 
pH environment, the three a-helices of the wild-type are unfolded as of the 12 14V and 
S173N mutants under neutral or low pH environment and at last should be unfolded into 
rich /3-sheet structures. 



Under low pH environment, the salt bridges of the wild-type (Table 1) are broken 
thus leads to the unfolding of the stable helical structures of RaPrP. The following salt 
bridges should contribute to the stabihty of RaPrP: ASP146-ARG147, GLU210-ARG207, 
GLU206-LYS203, GLU206-ARG207, ASP146-HIS139, GLU151-ARG150, GLU151- 
ARG147, GLU151-ARG155, ASP177-ARG163, GLU145-ARG135, ASP143-HIS139, 
GLU145-HIS139, GLU195-ARG155, ASP146-ARG150, ASP177-HIS176, ASP143-ARG147, 
GLU195-LYS193, HIS186-ARG155, HIS186-LYS184, ASP166-ARG163, GLU210-HIS176, 
HIS139-ARG135, GLU199-LYS193, GLU206-HIS176. Among these salt bridges, com- 
pared with the salt bridges of I214V mutant and S173N mutant listed in Table 1 and the 
observation from Figures 1~3, we find that GLU-195-ARG155 and ASP166-ARG163 
are important to contribute to the structural stability of wild-type RaPrP. 

Zhang (2011) found that there always exist salt bridges between ASP202-ARG156, 
ASP178-ARG164 in human and mouse prion proteins, between ASP201-ARG155, ASP177- 
ARG163 in RaPrP EOl . The author broke salt bridges of of human, mouse and rabbit 
prion proteins by doing MD simulations from neutral to low pH environment; conse- 
quently, the secondary structures of human and mouse prion proteins were not changed 
very much, while the stable helical structure of wild-type RaPrP collapsed [20 1. This is to 
say that the human and mouse prion protein structures are not affected by removing these 
salt bridges but the structure of RaPrP is affected very much by these salt bridges. There- 
fore, compared with human and mouse prion proteins, the salt bridges such as ASP201- 
ARG155, ASP177-ARG163 contribute greatly to the structural stability of RaPrP 



5 A Concluding Remark 



Like the controversy on "prion" theory, there is also a big controversy over "whether 
rabbits are resistant to prion infection or not?". This paper briefly reviewed the research 
results on rabbit prion protein and rabbit prions. This paper also did MD simulation 
studies on RaPrP and its I214V and S173N mutants with a little preliminary analysis of 
the role of salt bridges that plays in the structural stability of RaPrP. Further studies to 
understand the specific mechanism of RaPrP are clearly needed. 
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Figure 1: (Seedl) Secondary Structure graphs for rabbit prion protein wild-type, I214V 
mutant and S 17 3N mutant (from up to down) (x-axis: time (0-30 ns), y-axis: residue 
number (124-228); left column: neutral pH, right column: low pH. 
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Figure 2: (Seed!) Secondary Structure graphs for rabbit prion protein wild-type, I214V 
mutant and S173N mutant (from up to down) (x-axis: time (0-30 ns), y-axis: residue 
number (124-228); left column: neutral pH, right column: low pH. 
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Figure 3: (Seed3) Secondary Structure graphs for rabbit prion protein wild-type, I214V 
mutant and S173N mutant (from up to down) (x-axis: time (0-30 ns), y-axis: residue 
number (124-228); left column: neutral pH, right column: low pH. 
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Table 1: Salt bridges of the wild-type and the I214V and S173N mutants under 
neutral pH environment for the MD simulations with occupied rates for seedl, 



eedi, where ^ c 


lenotes the occu 


Die 


a rate is less tl 


lan 1U%. 


Salt bridges \ 


wild-type 


12 14V mutant 


S173N mutant 


ASP146-ARG147 


(100, 100, 100) 


(100, 100, 99.90) 


(100, 99.97, 100) 


GLU210-ARG207 


(99.15,97.58,98.43) 


(35.08, 98.08, *) 


(99.25, 55.00, 100) 


GLU206-LYS203 


(97.67, 92.85, 87.35) 


(29.07,41.00,57.78) 


(*, 32.52, 46.10) 


GLU206-ARG207 


(84.25, 84.80, 75.50) 


(65.78, 57.03, 73.68) 


(47.33,27.07,18.75) 


ASP146-HIS139 


(*, 73.65, *) 


(69.13,*,*) 


(*, 38.93, *) 


GLU151-ARG150 


(72.13,47.10,84.25) 


(65.35, 68.42, 88.25) 


(57.32,67.85,75.15) 


GLU151-ARG147 


(63.13,36.28,64.22) 


(39.02, 37.32, 34.45) 


(*,18.67, 73.33) 


GLU151-ARG155 


(60.78, *, 60.07) 




(*, *, 29.00) 


ASP177-ARG163 


(47.80,40.38,21.92) 


(*,25.17, *) 




GLU145-ARG135 


(33.52, *, 47.13) 


(*, *, 13.30) 




ASP143-HIS139 


(27.95, *, 10.95) 


(*, *, 29.40) 


(*, 26.97, 13.53) 


GLU145-HIS139 


(25.00, *, 18.60) 




(*, *, 26.72) 


GLU195-ARG155 


(21.88,*,*) 






ASP146-ARG150 


(21.23,66.38,34.25) 


(59.57, 15.38, 29.25) 


(42.77, 66.85, 12.42) 


ASP177-HIS176 


(21.17, 16.55, 10.15) 


(70.05,42.13,20.30) 


(25.53, 34.40, *) 


ASP143-ARG147 


(19.73, 50.30, 28.38) 


(33.85, *, 21.23) 


(10.43,42.20,61.53) 


GLU195-LYS193 


(19.68, 21.98, 10.27) 


(*, 19.53, *) 


(26.15,20.60,*) 


HIS186-ARG155 


(14.01, 96.35, 12.47) 


(*,*, 11.37) 




HIS186-LYS184 


(*, 37.95, 12.95) 


(*, 30.79, *) 


(82.28, 86.89, *) 


ASP166-ARG163 


(*, *, 54.08) 






GLU210-HIS176 


(*, *, 33.37) 




(24.25, *, *) 


HIS139-ARG135 


(*, *, 32.02) 




(*, *, 54.32) 


GLU199-LYS193 


(*, *, 11.95) 


(17.77, *, *) 


(*, 18.33, *) 


GLU206-HIS176 


(*, *, 11.83) 


(*, 62.37, *) 




HIS176-ARG163 




(72.11,*,*) 


(*, 39.57, *) 


ASP201-LYS203 




(32.53,87.10,91.17) 


(38.32, 69.48, 73.22) 


GLU199-LYS184 




(25.53, *, *) 


(18.22, *, *) 


GLU195-HIS186 




(23.32, *, *) 


(*, 25.82, 16.67) 


HIS186-LYS193 




(27.65, *, *) 


(*, 10.70, *) 


GLU145-ARG147 




(*, 80.58, 24.40) 




GLU210-HIS139 




(* 


28.42, *) 




ASP201-ARG155 




(* 


27.60, *) 




HIS176-HIS186 




(* 


20.18,*) 




GLU220-ARG227 




(* 


10.60, *) 




GLU206-LYS184 




(* 


*, 46.75) 




ASP177-ARG207 




(* 


*, 40.75) 




HIS176-ARG207 




(* 


*, 15.67) 




GLU151-LYS193 




(* 


*, 14.55) 




GLU206-HIS186 




(* 


*, 13.08) 




ASP143-ARG135 






(24.60, *, *) 


GLU199-HIS186 






(18.68, *, *) 


GLU206-ARG155 






(17.50, *, *) 


ASP201-HIS186 






(14.48, *, *) 


GLU145-ARG150 






(13.33, *, *) 


HIS176-LYS184 






(*, 21.93, *) 


GLU195-ARG227 






(*, 12.28, *) 


ASP146-ARG135 






(*, 10.70, *) 


HIS176-ARG227 






(*,*, 11.12) 


GLU220-HIS176 






(*,*, 24.17) 


ASP201-ARG207 






(*, *, 15.55) 


GLU210-ARG227 






(*, *, 12.00) 
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